Creation on demand of higher orbital states in a vibrating optical lattice 
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It is shown that the extended Hubbard Hamiltonian describing atoms confined in an optical lattice 
always contains commonly neglected terms which can significantly change the dynamical properties 
of the system. Particularly for bosonic systems, they can be exploited for creating orbital states 
on demand via the parametric resonance phenomenon. This indicates an additional application 
for optical lattices, namely the study and emulation of interactions between particles and lattice 
vibrations. 
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The past decade of experiments on cold atom systems 
confined in optical lattices has brought huge progress 
in simulating different variations of the Hubbard model 
P, 01 ■ It has become possible to experimentally mimic 
not only the simplest Bose-Hubbard model [sl, B], but 
also to study systems with internal degrees of freedom 
Recent theoretical works 0413 ^'^d first exper- 
iments considering higher bands in optical lattices |lll[l2| 
have opened a new, promising area of ultra cold atom 
research - orbital physics Typically, such systems 

are described with Hubbard-like Hamiltonians extended 
by additional inter-orbital interaction terms. In previous 
analyses, some of these additional terms, due to their non 
resonant character, have always been treated as unimpor- 
tant and therefore omitted. In this letter, I give a simple 
counter-example that shows that in the case of oscillat- 
ing optical lattices they can efficiently transfer atoms to 
liigher bands in a fully controlled way. All the argu- 
ments presented here are given directly for bosons inter- 
acting via 5-Vike interactions. Nevertheless, the central 
observation made here originates from the structure of 
any many-body Hamiltonian. Therefore, it can be easily 
adopted for other interparticle interactions as well as for 
fermions. 



For simplicity, let me consider spinless (or polar- 
ized) bosons interacting via two body (5-like contact in- 
teractions confined in a 2D optical lattice T4xt(^) — 
QxSvci'ikx) + qySvcL'{ky) -t- ^^^z^, where k — 211 jX 
is the wave vector of the laser field and a;^ is the 
frequency of the trapping harmonic potential in the 
z direction. The optical lattices need not have the 
same depth in the x and y directions and they can 
be changed independently. The general Hamiltonian of 
this system written in second quantization formalism 

has the form -H = li?r^\r) [-^^V^ -|- Vext(r-)j *(r)-l- 

£ jd^r 5'^(r)vl/t(r)xp(r)vl/(r), where ^(t-) is a bosonic 
field operator and g is the strength of the contact in- 
teractions. It is convenient to measure lengths in units 
of the laser wavelength A, and all energies in units of 

the recoil energy — ^2ni\'^ • '^^^ trapping poten- 
tial in the z direction is characterized by the dimen- 



sionless quantity k — TiuJz/'^Eji. The dimensionless cou- 
pling constant is g = 167r^ao/A, where ao is tire s-wave 
scattering length. In typical experiments with ®^Rb or 
^^Cr confined in an optical lattice, far from the Fesh- 
bach resonance, 17 ~ 1. To derive an extended Bose- 
Hubbard model describing this system, one expands the 
field operator in the ground and first excited Bloch bands 
as follows: \l/(r) "i'/'i + bi4>f{r) + Ci4if{r). 

The functions (j)°{r) = {x)y° (y) Z (z) , 0f(r) = 
Xl{x)y^{y)Z{z), and 0f(r) = A:j^{x)yl{y)Z{z) are 
products of one-dimensional Wannier functions A'"(x) 
(3^f (2/)) from band a localized in site i, and the ground 
state of the harmonic oscillator in the 2 direction Z{z) — 
(K/7r)^/^exp(— K2^/2). This decomposition is valid pro- 
vided that ^ Qx/Eii and ^ qy/En. In such a 
case, the energy gap for excitations in the z direction is 
much larger than the gaps in the lattice directions, and 
therefore particle dynamics is frozen in the z direction. 
Bosonic operators a^, 6^, and Ci annihilate particles at 
site i in the s, p^-, andpj, orbitals respectively. For conve- 
nience, I introduce particle number operators ni*^ = alai, 
fix '' — blbi, fly '' — cjci as well as the dimensionless alge- 
braic vector Q = [q^/ Efi,qy / Ee.t n) which characterizes 
tlie geometry of the optical lattice. In the above approx- 
imation, the Bose-Hubbard Hamiltonian takes the form: 



- E J(ii4<^, + 4bj) + jfcic„ (lE 



where is the standard one dimensional nearest neigh- 
bor hopping amplitude in the direction d for band a. The 
summation a is understood as a summation over 

all sites i and over all nearest neighbors j of site i in the 
direction d. The on-site Hamiltonian Hi is a sum of the 
single particle energies and two-body interaction terms 
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FIG. 1: Parameters of the Hamiltonian ([T]) as functions of 
lattice depth for the symmetric case Q — (q, q, 8) and g = 1. 
Parameters originating in contact interactions are at least ten 
times smaller than the energy gap between the s and p-bands. 



and has the form: 



E,hf + ^nUnf - I) 



J2 U^a'fifnf 
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(lb) 



These summations run over orbital index a € {s,x,y}. 
The single particle energies E^- depend only on the lattice 
geometry while all the parameters U depend additionally 
on the dimensionless coupling g. They can be calculated 
directly: 
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V2 + V;xt(r) 



(2) 



To get more insight into the system, the values of 
these parameters for the example of a symmetric lattice 
Q = ill 9; 8) and g = 1 are presented in Fig[TJ For exper- 
imentally available systems, far from the Feschbach reso- 
nance, all contact energies U are at least ten times smaller 
than single-particle excitation energies AEg- = E^ — Eg 
where a € {x,y}. Usually, one can neglect the two next 
to last terms in the Hamiltonian (jlbp , if the contact in- 
teraction energy is a small correction to the gap energy 
between bands. These terms describe interaction pro- 
cesses which transfer two particles between bands and 
they are highly suppressed since they violate energy con- 
servation. However, if one considers a scenario when the 
lattice parameters Q or the coupling constant g vary in 
time, that argument is no longer valid since energy con- 
servation simply does not hold in such a case. Never- 
theless, in previous analyses 0, S these non resonant 
terms were always neglected, even for fast- varying Hamil- 
tonians (l3 - [l6| . Under such a far-fetched approximation, 
particle numbers in each orbital are constant. In this pa- 
per we utilize these commonly neglected terms and pro- 
pose a mechanism for creating orbital states on demand. 
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FIG. 2: Creation of the p-band states by a vibrating optical 
lattice with Q = (32 + 4sm{ujt), 20, 8). Plot (a) presents the 
transfer efficiency as a function of vibration frequency uj. Two 
distinct peaks are visible. For these particular frequencies, p- 
band states become highly occupied. The dashed line comes 
from a generalized model that also takes into account d-band 
orbitals. Additional peaks correspond to resonant frequencies 
in which one of interacting atoms can be promoted to the dx 
or dy band, respectively. Moreover, resonant frequencies cji 
and ijj2 are shifted. Plot (b) presents the occupation of basis 
states as a function of time for the resonant frequencies. At 
frequency uji {iJ2) the py (px) orbital is populated. Plot (c) 
shows a comparison between the situation when tunneling is 
totally neglected (solid line) and when it is taken into account 
(filled circles). The inset shows the variance of the on-site 
number operator as a function of time for cj = a;2. Detailed 
explanations are given in the text. 



To show that this scenario is realistic in present day 
experiments, numerical simulations are shown for ^^Cr 
atoms confined in an optical lattice with A = 523 nm. 
The contact interaction coupling is g « 1.8. 

The starting point of the analysis is the many-body 
ground state of bosons confined in a static and deep 
optical lattice Qq = (32, 20, 8) filled with two bosons 
per site on average. In such a case, tunneling is highly 
dominated by the on-site contact interactions, and the 
many body ground state is the Mott Insulator in the 
s-band. Hence the state of the system can be charac- 
terized quite well by a product of local ground states in 
independent lattice sites. Now let us study the situation 
when the lattice parameters change periodically in time 
in the following way: Q{t) — Qo + SQ{t). In the simplest 
case, the lattice parameters oscillate with some ampli- 
tude A and frequency uj in one spatial direction only, i.e. 
SQ{t) = {A, 0, 0) sin(a;i). Since the lattice depth is large 
for our choice of Q the tunneling processes are very slow. 
Therefore it will be assumed that the dynamics in each 
site is independently governed by the single-site Hamil- 
tonian Hi and the number of particles in each site is con- 
served. The correctness of this assumption was verified 
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with dynamical many-body calculations and is discussed 
below. At the initial moment, two particles occupy the 
s-band state. Therefore, due to the structure of the 
Hamiltonian (jlbp , the entire dynamics takes place in the 
subspace spanned by three states: |200) — ^a'^^jvac), 

|020) = ^fcl'^lvac), and |002) = ^ct2|vac). In this sub- 
space, the Hamiltonian has a simple matrix form 

l2Es + Uss Usx Usy \ 

H{Q{t)) = 2E, + Uxy . (3) 

\ Uxy '^-^y ~^ ^yy / 

All the parameters of this matrix depend on time through 
the time dependence of the lattice shape Q{t). To quan- 
tify the influence of the vibrating lattice on the state of 
the system we define the transfer efficiency as the high- 
est depletion of the initial state for a given frequency a; 
and amplitude A. In Fig. [2^ (solid line) this transfer 
efficiency is presented as a function of frequency cj for 
amplitude A — 4. It is clear that for two characteris- 
tic frequencies lji and uj2, the initial state can be totally 
depleted. Fig. [2}d presents the time dependence of oc- 
cupations for the corresponding two frequencies. The 
full width at half maximum for both resonances is about 
Slo /{2Tr) « 700 Hz. The characteristic frequencies almost 
do not depend on the amplitude A and they are approx- 
imately equal to the energy difference between the ap- 
propriate eigenstates of the matrix H{Qq). Let me note 
that a full transfer of interacting atoms is obtained in a 
few milliseconds and therefore it is much faster than the 
experimentally obtained decay time of hundreds of ms 

To show that the predictions described are almost in- 
sensitive to the approximations in the model, two ad- 
ditional tests were performed. Firstly, the expansion of 
the field operator '^{r) were generalized so that all d- 
orbital states and their possible interactions were taken 
into account. The resulting transfer efficiency is shown 
with a dashed line in Fig. [5^. As is seen, the previously 
predicted frequencies are slightly shifted and two addi- 
tional peaks have appeared. They correspond to the res- 
onant frequencies in which one of the interacting atoms 
is promoted to the dx or dj,-bands, respectively. This 
result shows that the described mechanism for creating 
higher orbital states is highly selective, and that for a 
particular choice of the resonant frequency one can ne- 
glect the other orbital states, i.e. the whole dynamics 
effectively takes place in a two dimensional subspace of 
coupled states. Secondly, to check the influence of the 
tunneling processes, the full many-body dynamics in a 
ID optical lattice was studied. Due to computational 
complexity, calculations were performed for eight atoms 
in a lattice with four sites, and with periodic boundary 
conditions. Vibrations Q = (32 -I- 4 sin(a;t), 20, 8) with a 
frequency close to W2 were considered. The simulations 
performed evidently show that for the chosen lattice pa- 
rameters, the dynamics of the system can be treated as a 
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FIG. 3: Two experimental scenarios for the creation of a two 
particle superposition in excited band states. Initially the 
system is prepared in the ground state of a highly non sym- 
metric lattice Qo = (32,20,8). Plot (a): In the first sce- 
nario, two particles are transfered to the Py state by apply- 
ing appropriate vibrations (Pi) and then lattice depths in 
both directions are equilibrated (E). If this process is slow 
enough, then the final state of the system is in the superpo- 
sition (1 020) - 1 002)) 7^2. Plot (b): In the second case, two 
particles are transfered to the py state by applying appropri- 
ate vibrations (Pi). When the initial state is half-depleted 
the vibration frequency is changed to the other resonance, 
and the px state is filled (P2). Then similarly to the previous 
scenario, lattice depths are equilibrated (E). The final state is 
in a complex superposition of basis states and the occupation 
of each basis state varies in time (F). At the moments when 
the occupation of px and py orbitals are equal, the system is 
in one of the vortex states (|020) ± il002))/^. 



dynamics carried out independently in each lattice site. 
Fig. [2J3 shows the resulting transfer efficiency (filled cir- 
cles) compared to single-site predictions (solid line) as 
well as the variance of the on-site number operator for 
the resonant frequency UJ2 (inset). The situation changes 
significantly for larger tunneling amplitudes, i.e. for shal- 
low lattices (about I2i?/{ for the example studied), and 
it will be discussed elsewhere. 

Additionally, let me briefly discuss the appreciably 
more complicated situation when the initial lattice is 
symmetric in both directions. Then, both p-band basis 
states have the same energy, but due to the existence of 
the coupling Uxy they are not eigenstates of the Hamilto- 
nian. As previously, resonant frequencies are determined 
by the eigenenergies of the Hamiltonian ([3]) but now par- 
ticles can be excited independently to symmetric or an- 
tisymmetric combinations of the original states |±) — 
(|020) ± |002))/a/2 by applying symmetric or antisym- 
metric lattice vibrations 6Qj.{t) = {A, ztA, 0) sin(wi), re- 
spectively. 

The mechanism discussed above can be easily extended 
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to the idea of creating an orbital superposition of states 
(1 020) + e''^|002))/\/2 with an arbitrarily chosen phase 
(j). This can be done in a variety of different ways. Let 
me consider two scenarios which are efficient for creat- 
ing states with = kTT/2, where k — 0, . . . , 3 (see Fig. 
[3]). In both scenarios one starts with the system pre- 
pared in the static, non symmetric lattice with parame- 
ters Qq = (32, 20, 8) in the insulating ground state with 
two particles in each lattice site. Then, the system is 
driven by a vibrating lattice with frequency wi (^2) and 
the state |002) (|020)) becomes occupied. In the first sce- 
nario (Fig. [3^) one completly depletes the ground state 
(time interval Pi). In the second one (Fig. [Sb) the driv- 
ing frequency is switched to W2 when the initial state has 
become half-depleted. Vibrations with frequency uj2 lead 
to the occupation of the state |002) (P2)- In both sce- 
narios, when the ground state becomes totally depleted, 
the lattice depth is brought down to equate Qx and qy 
(interval E). From this moment, the basis states |020) 
and 1 002) have the same energy, but due to the contact 
interactions (the last term in the Hamiltonian (jlbp ) they 
are not eigenstates of the Hamiltonian. The final state 
of the system (F) is in some superposition of the basis 
states. The relative phase depends on the speed and 
details of the equilibrating process. Nevertheless, in the 
extreme case of adiabatic equilibration, the phase rela- 
tionship between the orbital states is known. 

In the first scenario at the beginning of the equilibra- 
tion interval E the system is in an eigenstate (a ground or 
excited state in p-orbital subspace) of the initial Hamil- 
tonian ([3]). Therefore, during an adiabatic equilibration 
of lattice parameters the system has to remain in that 
eigenstate of the temporal Hamiltonian. At the final mo- 
ment, the states are ^(|020)± |002)) (the relative phase 
is or tt). a similar situation occurs in the second sce- 
nario. Before equilibration the system is in an almost 
equal superposition of states |020) and e**|002), which 
are the eigenstates of the Hamiltonian at that time. The 
phase $ is very hard to control and therefore in practice 
it will be unknown. During the adiabatic equilibration 
the eigenstates track the time-evolving Hamiltonian, and 
change to the pair (|020) ± |002))/y2 in terms of basis 
states of the final Hamiltonian. During this process the 
eigenstates accumulate phases A+ and A_, respectively. 
It can be shown straightforwardly that the final state 
of the system is e'" [cos(/3)|020) -I- i sin(/3)|002)], where 
a = {A++A_ + $)/2 and /3 = A+ - A_ - $. At later 
times (interval F) the relative phase factor between |020) 
and 1 002) remains zLi, while [3{t) evolves and leads to 
Rabi oscillations. At the moments when occupations of 
basis states are equal the state is in one of the vortex 
states (|020) ± i\002))/y/2 (see Fig. Oj). 

Let me note that in the model discussed here, exita- 
tions to higher orbital states take place in all sites of the 
optical lattice not only independently but also simulta- 
neously. Nevertheless, by applying an additional weak. 



anharmonic external potential one can make the energy 
gap between the orbitals site dependent, and, in conse- 
quence, the resonant frequencies can differ from site to 
site. This gives not only temporal but also spatial con- 
trol of excitations. As such, it could be utilized for quan- 
tum computation engineering to address selected qubits 
formed by lattice sites with two bosons 

The results presented here show that lattice vibrations 
can effectively couple the ground Wannier state to cho- 
sen orbital states. The mechanism is very general since 
it originates from the fundamental Hubbard Hamiltonian 
by taking into account interactions between particles. In 
general, for each term describing a non vanishing interac- 
tion energy between particles in different single-particle 
states there always exist corresponding term describ- 
ing a transfer of particles between them. In the static 
case these terms usually violate energy conservation and 
can therefore be neglected. However, when one consid- 
ers time dependent Hamiltonians they should be always 
taken into account because they can significantly change 
the dynamics of the system. The additional terms in the 
Hamiltonian that were studied here should also be taken 
into account when rapid changes of the Hamiltonian are 
considered, e.g. to correctly describe quantum quench 
problems 

Finally, let me note that instead of driving the sys- 
tem via a changing depth of the lattice one can alter- 
natively change the aspect ratio k or coupling constant 
g. The last possibility is quite interesting since it can 
be realized by applying an oscillating external magnetic 
field in the proximity of a Feshbach resonance. It means 
that by preparing appropriate pulses of the external mag- 
netic field one could mimic interactions between particles 
and phonons by artificial vibrations propagating along 
the lattice (artificial phonons). Thus, the mechanism de- 
scribed here is not only an additional way of playing with 
orbital physics but it may also lead to a better under- 
standing of those solid state problems that are beyond 
static theory [l9| . 
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